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fSJ ■ Abstract 

o: 

r^ . Fluctuations of the mean field of a globally coupled dynamical systems 

t^ I are discussed. The origin of hidden coherence is related with the instabil- 

_!. ■ ity of the fixed point solution of the self-consistent Perron-Frobenius equa- 

OS ■ tion. Collective dynamics in globally coupled tent maps are re-examined, 

pj I both with the help of direct simulation and the Perron-Frobenius equa- 

^^' tion. Collective chaos in a single band state, and bifurcation against initial 

conditions in a two-band state are clarified with the return maps of the 

O I mean-field, Lyapunov spectra, and also the newly introduced Lyapunov 

(— I ' exponent for the Perron-Frobenius equation. Future problems on the col- 

O ■ lective dynamics are discussed. 
>'■ 
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Between the first and last workshops at Como, we have seen some progress in the 
study of globally coupled chaotic systems. In globally coupled maps, the studies 
have elucidated the clustering0, the breakdown of the law of large numbers 0, 
chaotic itinerancy P, H, ^, complexity in the partition to clusters 0, the exis- 
tence of high-dimensional tori |^, ^, 18|, and marginal stability with information 



avalanches [pTq] . The universality of these behaviors was confirmed by a variety 



of mappings 0, H, |^, as well as ordinary differential equations [|T0|- Related 
(but partially different) behaviors have been studied in oscillator-type systems 
J^, |, |1^, [1^, |I9|, globally coupled Ginzburg-Landau equations |]1T], [T2|, 0, and 
globally coupled tent maps 0, [l^ (see also [jl6[). Applications to Josephson 
junction arrays []^, |TB], ^ and optical systems have also been disussedpU]], while 
relevance of the above ideas to biological networks has been studied together with 
extensions to inhomogeneous systems [^, |2l|, |2^, |23l . 

It is noted that synchronization of globally coupled oscillators has been studied 



since the 70's [Ba, Eq, 27, 29, 28|. However, there is an essential difference between 



these earlier studies and recent studies since 1989. The concern of the former was 
synchronization in a system with (randomly) distributed frequencies (or other pa- 
rameters) without explicit interest in complex dynamics (such as chaos), while 
the main focus of the latter has been the interplay between chaos and synchro- 
nization, where identical elements can be differentiated and desynchronized due 
to the chaotic instability. 

The simplest model among such globally coupled dynamical systems is the 
globally coupled map (GCM), originally introduced (0 see also 0) as a mean- 
field- type extension of the coupled map lattices (CML) [^. The GCM is a 
discrete-time dynamical system whose elements interact with all other elements. 

Here we focus on the following class of models: 

e ^ 
x.+i(z) = (1 - e)/(x„(z)) + -Y: /(x„(j)), (1) 

^^ i=i 
where n is the discrete time and i is the index of elements (i = 1,2, ■ ■ ■ , A^=system 
size)P]. The mapping function /(x) is chosen so that it can show chaos. In the 
present paper we mainly study the logistic map 

fix) = 1 - ax^, (2) 

or the tent map 

f{x)=a{0.5-\x-0.5\), (3) 

since there have been investigated in detail as prototypes of dissipative chaos. 

Earlier we reported the breakdown of the law of large numbers and the emer- 
gence of subtle coherence [^ in the desynchronized (turbulent) state. Here the 
temporal fluctuations of the mean-field 

/^n = ^E/(^W) (4) 

i 

are measured with the increase of the system size. The fluctuations turn out 
to decrease only up to some size, beyond which they remain constant. This is in 
contrast with the fluctuations in an ensemble of decorrelated random elements, 
where the mean square deviation of the mean- field h should decrease with A^~^. 
The above remnance of the mean field fluctuations, thus, implies a hidden co- 
herence among the oscillations of the elements. It should be noted here that 
this coherence is subtle in nature, and is neither due to (i) collective behaviors, 
nor (ii) the frozen relationships among elements. In an earlier paper we have 
attributed the origin of the breakdown of the law of large numbers to the insta- 
bility of the fixed point solution of the Perron-Frobenius (PF) operator. Here we 
discuss the relationship in a little more detail (see also ||3^ and ||l^). 



Before going into detailed discussions, we note that the interaction term in 
(1) could be replaced by noise (whose root mean square is 0(l/viV)), if the 
law of large numbers ( and the central limit theorem) were valid. Then, in the 
limit A^ ^ oo, the dynamics of (1) would be reduced to A^-independent logistic 
maps given by Xn+i{i) = (1 — e)f{xn{i)) + eh* , where h* = (l/N) J2i fi^ii)) is a 
constant mean field. 



2 Instability of the fixed point of the Perron- 
Frobenius Operator 

In the limit of A^ — > oo, the snapshot distribution Pn{y) is defined by the prob- 
abihty that x{i) assumes the value y {pn{y) = {^/N)Ylf=iS{y — x„(j))). This 



probability distribution can be calculated with the use of PF operator [|^, ^ 
It is written as [| 

where the two preimages are given by 

ypre = f ( YZ ''' ^^^ 



and f~^{y) = ±J^^ for our logistic map. As was discussed in |^, the breakdown 
of the law of large numbers may be associated with the instability of the fixed 
point solution of the self-consistent Perron- Frobenius (PF) equation. Indeed the 
instability can be explained as follows: 

If the law of large numbers is valid, the mean- field is constant, and thus the 
above distribution is time-independent, and should be given by the invariant 
measure p*{y) for the logistic map 

a;„+i(0 = (l-e)/(x„(z)) + e/i*, (7) 

where h* is determined by the self-consistent condition h* = J f{x)p*{x)dx. Since 
equation (7) is just a one-dimensional map, the invariant measure is straightfor- 
wardly calculated with the help of the Perron- Frobenius operator for the one- 
dimensional map. As is discussed in [H, the breakdown of the law of large num- 
bers implies that the above self-consistent probability solution either does not 
exist at all, or is unstable in case it does exist. Accordingly we have to study the 
existence and stability of the fixed point solution p*{x). 
For this study we first consider the one-dimensional map 



''ini 



a;n+i(0 = (1 - ^)f{xn{i)) + thi 

where hin is an "input mean-field" given as a parameter (which should be deter- 
mined later self-consistently). Then we compute the invariant measure p\^ {x) 
numerically, by using the Perron- Frobenius operator for the one-dimensional map. 
From this p^. {x) we obtain the "output" mean- field 



hout = j f{x)pl^^{x)dx. (9) 

The fixed point solution of our original self-consistent equation of course 
should obey /ij„ = hout- The obtained invariant distribution p*{x) from the 
one-dimensional map surely satisfies the self-consistent equation. 

^There is a typo in eq.4 of 0], which should be read as eq.(5) here. 



Instead of investigating the stability of p*{x) with regard to the original PF 
equation, we check the stability of the self-consistent solution against a perturba- 
tion of the mean field value only. If this solution is stable, any slight perturbation 
to p*{x) should decay with time. Aside from an exceptional form of perturbation, 
the value of the mean-field h is altered by it. Thus, the stability of the fixed point 
solution requires \dhout/dhin\ < 1 at /ij„ = hout- In Fig-1, we have plotted hout 
as a function of hm- As can be seen in Fig. lb), the condition \dhout/dhin\ < 1 is 
not satisfied here, which explains the origin of the breakdown of the law of large 
numbers. 

The complicated behavior in Fig.l is due to the existence of windows in the 
logistic map. Since windows exist in any neighborhood of the logistic parameter 
leading to chaos, the complicated behavior is expected to be a general feature. 
Thus we expect that the law of large numbers is not satisfied as long as the local 
mapping has windows. 

Fig.l 



In the tent map, on the other hand, the relationship between hout and /ij„ is 
smooth, and \dhoutldhin\ < 1 at the unique fixed point /ij„ = hout (see Fig. 2). 
This confirms the stability of the fixed point solution of the PF equation for 
the tent case, while a rigorous proof needs further studies. Indeed it has been 
numerically shown that the law of large numbers is valid in this parameter regime 
of globally coupled tent maps @ ( we call this regime completely desynchronized 
state). 

Fig.2 



Although the Perron-Frobenius approach may be powerful, the following ques- 
tions remain to be solved; 

• How can the crossover size Nc and the remaining correlations be determined 
from the PF equation? 

Intuitively, the crossover size may be related with the degree of instability 
of the fixed point solution. So far we have no formulation justifying this 
intuition, based on the PF equation. Numerically it has been found that 
the variance < {5hY >oo (at the infinity size limit) is inversely proportional 
to 1/A^c, and that it is proportional to e^ for small coupling e 0. 

• Dependence on noise 

By adding a noise larger than some threshold, the law of large numbers is 
known to be recovered. The threshold seems to be roughly proportional to 
< {5hY >oo- The probability distribution p in the stochastic case is written 

as 



Pn-liVp 



,„(.). _/..P„^E^_^^ ,0) 



with 



_ ,-u X - w -ej pn-i{z)f{z)dz ^ 

ypre,w J { T J) \^^) 

1 — e 

by using P{w) as the probabihty of the noise to assume the value w. Indeed 
Pn{x) approaches a fixed point distribution numerically, when the noise is 
large enough, as is shown in Fig. 3. This agreement is a starting point for 
further studies on the effect of noise in terms of the PF equation. 

Another mystery remaining to be solved is the anomalous recovery of the 
law of large numbers: When the noise amplitude is not large enough, the 
variance of the mean field decays with 1/N^ with /3 < 1 0. This anomalous 
power may be due to hierarchical destruction of coherence, which has to be 
unveiled in the future. 

Another striking mean-field behavior is seen for a GCM with non-identical, 
distributed parameters 0. Here the mean field fiuctuation decays with 
1/yN up to some size, but then increases to reach the original level of 
the non-distributed parameter, with the further increase of size. So far 
the origin is totally unexplained, although it is natural to expect that the 
windows existing everywhere in the parameter space at least play some role. 



-Fig.3 



3 Single-band State: High-dimensional Collec- 
tive Chaos 

In the globally coupled logistic map, transitions among coherent, ordered (clus- 
tered), partially ordered, and turbulent (desynchronized) phases are observed, 
which is rather common in globally coupled dynamical systems. In the globally 
coupled tent map, the coherent state is stable if (1 — e)a < 1 is satisfied, while 
we have not observed clustered states. 

Instead, we have observed desynchronized states which keep some coherence. 
Earlier we reported that the desynchronized states consist of single band, two- 
band, and completely desynchronized states |Q. 

In a (single) tent map an orbit is always expanding in contrast with the logistic 
map. The motion is linear in each branch ( here we call x < 0.5 the left branch 
and X > 0.5 the right branch). Although there are no windows in the chaotic 
region, there are band splittings. For 1 < a < \pl, Xn takes a value larger or 



smaller than x* = a/(l + a) (the unstable fixed point), alternately in time. Hence 
two orbits starting from x > x* and x < x* do not mix, and the Perron-Frobenius 
equation has a period-2 solution there. 

Thus the stability argument of the fixed point in the previous section does not 
hold any more, if the self-consistent tent map Xn+i = (1— e)a(0.5— |x— 0.5|)+e/i* = 
F{x) is in a two-band region (where h* is determined self-consistently as in the 
previous section). The map F{x) must be in a single-band region, to show a 
completely desynchronized state, i.e., a and e should satisfy the condition given 
by F{F{F{l/2))) > x*. For smaller a or larger e where the condition is not 
satisfied, the map shows two-band splitting, where the completely desynchronized 
state disappears and the globally coupled map yields the two-band motion as 
mentioned in [Q]. 

In the present section, we study the mean field dynamics of the single-band 
state in more detail. In the single-band state, all elements oscillate almost coher- 
ently, that is, elements always fall within a narrow range (see Fig. 18 (b) of 0). 
Thus the snapshot distribution p„(x) is finite only within a narrow interval of 
x, and the amplitude of the mean field dynamics does not decay at all with the 
increase of the system size. 

Here, it is rather straightforward to expect some collective motion, represented 
by the mean-field dynamics. In Fig.4, the return map of the mean field is plotted, 
obtained from direct simulation of a GCM with A^ = 50000. We note that the 
return map from the PF equation agrees well with Fig.4, and the scattered points 
around x = 0.5 are not due to a finite size effect. For x > 0.56 and x < 0.45, 
the plot exactly agrees with f{x). Indeed, if all the elements fall in the same 
branch (according to x >< 0.5), J2j fi^U)) = fiJ^j ^U))y since f{x) is linear for 
each branch. Thus we get hn+i = f{hn), when all elements fall within the same 
segment of the tent map ( This statement holds generally for a globally coupled 
piecewise-linear map). Since the elements are not in complete synchronization 
(although close), they split into two branches when x„(i) is close to 0.5, leading 
to the deviation of the return map from f{x) around 0.5, shown in Fig.4. Still 
we should note that the mean-field follows the same dynamics as in the coherent 
case for most of the time steps n with Xn distant from 0.5. The "almost" coherent 
behavior arises from this use of the same branch. 

Fig.4 



Correlated motion of elements within a small band leads to collective dynam- 
ics. From the time series of the mean field, we have measured the correlation 
dimension by increasing the embedding dimension, expecting to detect a low- 
dimensional attractor. So far, however, a saturation in the slope of the corre- 
lation integral P^ is found neither from the direct simulation of the GCM nor 



from the PF equation. Thus the dynamics of the mean field is still very high. 

Lyapunov spectra, on the other hand, have a distinct character. As is shown 
in Fig. 5, one exponent is much larger than the other Lyapunov exponents, all of 
which are positive. Furthermore, the number of distinct exponents remains one, 

6 



even if the number of the elements is increased. 

Usually, in spatially extended systems, Lyapunov spectra have a scaling behav- 



ior such that Aj = A{i/N) |^, 36|. This is also true in the desynchronized state 
of the globally coupled logistic map . Recently Nakagawa and Kuramoto have 
found that in a desynchronized state of the globally coupled complex Ginzburg 
Landau equation, only a few distinctly larger (positive) exponents exist, whose 
number remains constant with the increase of the system size. Other scaled 



exponents are null or slightly positive, depending on the parameter [12|. This 
characteristic feature of the Lyapunov spectra is similar to our findings for the 
single-band state here, although the remaining exponents are always positive in 
our case. 

Fig. 5 



4 Two-band State 

When the nonlinearity is further increased (or the coupling is decreased) in the 
globally coupled tent map, we have a two-band attractor besides the single-band 
attractor, as is discussed in [||. Here elements split into two groups, which 
oscillate out of phases with each other. If the oscillation of the elements within 
the same group were completely synchronized, this state would be a two-cluster 
attractor as studied in the globally coupled logistic map . Here elements within 
each group are not synchronized, although they oscillate within a narrow band, 
using the same branch. This separation is stable: Once the groups are formed, 
the constituents of each group are invariant. Thus, there are many attractors 
with different partitions, coded by the number of elements in each group ( for 
example 45:55, 40:60, 50:50 and so on, for A^ = 100) |. 

Now let us study the mean-field dynamics for this 2-band state. In Fig. 6, we 
have plotted return maps of the mean-field by changing parameters. We clearly 
see collective motions close to quasiperiodic behavior. 

The mean field dynamics changes its nature according to the number of ele- 
ments in each group. "Bifurcations" with the change of the partition is given in 
Fig. 7, where the return maps of 5 attractors with different partitions are plotted. 
As the imbalance between two groups is larger, the amplitude of the mean-field 
dynamics increases, and the motion shows a quasiperiodic behavior more clearly. 

We have also numerically solved the PF equation starting from a variety of 
initial distributions po{x). Again, depending on po{x), we have obtained a variety 
of spatial return maps corresponding to those obtained from the direct simulation 
(see Fig 8 for 3 examples, although patterns corresponding to other cases are also 
found). Note that the width of scattered points remains finite even in the PF 
equation ( i.e., the N ^ oo limit), and thus that the attractor is higher than 

^ Note again that clustering in the original sense is not seen in the globally coupled tent 
map. 



one-dimensional ( i.e., the section of a 2-diniensional torus). 

It should be noted that the dynamical state "bifurcates" with the change of 
the number of elements. This "bifurcation" is quite analogous to the cluster 
bifurcations found in most GCM [^ . Indeed, if all the elements within the same 
group were synchronized, our model would reduce to the two dimensional map 

XUi = (1 - ^)f{K) +eJ:iN,/N)fiXi). (12) 

Then the bifurcation with Nj would be nothing but that with the change of 
the coupling parameter in the two-dimensional map. Since elements in the same 
group are not synchronized but agree only within a band width 6, the above 
reduction is possible only within the precision 6, and the cluster bifurcation serves 
only as an approximate explanation of the change of the return map. 

The Lyapunov spectra do not show any remarkable change against the "bi- 
furcation" with the number of elements in a group. All exponents take rather 
close values, as shown in Fig.5 a) and Fig 9a). Contrary to the clustering case 
0, a stepwise structure is not observed either, which might be expected if the 
oscillation in each group were synchronized ( see Fig. 9). Thus the above ap- 
proximation to two clusters is rather crude, and the deviation in each group is 
not negligible for the total dynamics. We also note that, when the single-band 
and 2-band states coexist, the former has a slightly larger maximal Lyapunov 
exponent (see Fig. 9b), while the Kolmogorov-Sinai entropy does not show any 
difference between the two states, within our numerical accuracy. 

Although collective motion can clearly be seen, we have not yet found a case 
where the mean-field dynamics clearly shows a low-dimensional torus. The slope 
of the correlation dimension increases beyond 4, and it is hard to confirm the 
low-dimensional motion. 

-Fig.6 

Fig. 7 



-Fig. 8 
-Fig. 9 



5 Lyapunov exponent for the Perron-Frobenius 
equation 

In order to investigate the stability of the mean-field dynamics, we consider the 
evolution in a tangential space of the probability function p„(x). First apply a 
tiny perturbation Sp{x) to make p'n{x) = Pn{x)+Spn{x) (/ 6p{x) = to satisfy the 
normalization condition) and see how this 5pn{x) evolves. Applying a standard 
method for calculating the Lyapunov exponent, we measure the growth rate by^ 



hmr^^hms,,^o E f^og^rj^-^^^^. (13) 

In our self-consistent PF equation, the above exponent is negative, when pn{x) 
approaches a fixed point (and the mean field satisfies the law of large numbers). 
On the other hand, the exponent is positive, when the collective dynamics is 
chaotic. 

The change of the Lyapunov exponents with parameters a and e is given in Fig. 
10. With the increase of e, bifurcations from the completely desynchronized state, 
and the 2-band state with collective torus, to the (high-dimensional) collective 
chaos is seen as the change of the sign of the Lyapunov exponent. From our 
numerical results by fixing e at .3, the dependence of the Lyapunov exponent 
is roughly proportional to (oc — a) (for a < ac), with ac ~ 1.81, beyond which 
the quasiperiodic motion of mean-field dynamics is seen. Here the existence of 
another branch ( for a < 1.6 ) is due to the single-band motion. It is shown that 
the single band motion gives a much smaller Lyapunov exponent. In principle it 
must be possible to detect some variation against the partition to groups in the 
2-band state, which, so far, has not clearly been confirmed, as in the case of the 
Lyapunov spectra from the GCM in §4 Q. 

Fig. 10 

The Lyapunov spectra of the PF equation may be a useful tool to analyze the 
instability of the mean-field. We have several problems to be studied in future: 

• Here we have discussed only the maximal exponent, but measuring all the 
exponents will be important in order to detect the high- dimensional collec- 
tive motions, which are otherwise hard to confirm. 

• Possible relationships between the Lyapunov spectra in the original GCM 
and the spectra in MF: The relationship, which one may expect to exist, 
cannot be so straightforward. For example, in the completely desynchro- 
nized state, all the exponents in the GCM are positive, but the maximal 
exponent in the PF equation is negative. When the collective motion shows 
quasiperiodic motion with a null exponent in the PF, the exponents in GCM 
are again positive. Nontrivial relationships between two Lyapunov spectra 
should be explored, for which detailed studies of A^-dependence of the spec- 
tra in GCM may be necessary. 

• Since the exponent from PF equation gives the (in)stability of the mean- 
field dynamics, one may expect to find some answers for the remaining 



■^For computation, J Spn{x)^dx is rescaled to a given tiny value ( say 10 —5) at every time 
step like in the usual computation of Lyapunov exponents without the use of the Jacobi matrix. 

^ In Fig. 10, only the exponents from 5 initial distributions are measured. By sampling more 
initial distributions, we may expect changes of Lyapunov exponents according to the partition 
into two groups in the 2-band state. 



questions in §2. 

• Of course, it is interesting to describe the bifurcations of the mean-field 
dynamics ( against a and e as well as the partition to groups) in terms of 
the Lyapunov exponents of the PF equation. 

6 Discussion 

In the present paper we have discussed the mean field dynamics of globally cou- 
pled maps. Some of the earlier reports are discussed with the use of the Perron- 
Frobenius equation. The mean-field dynamics in globally coupled tent maps 
is studied in more detail. The complete desynchronization in one phase is ex- 
plained by the stability of the fixed point of the PF equation, while it shows 
high-dimensional collective chaos in a single band state. In the two-band state 
the bifurcation from high- dimensional tori to chaos is observed with respect to 
the number of elements in grouping. Lyapunov analyses from the direct GCM 
and from the PF equation are presented. 

The globally coupled tent map has some features which distinguish it from the 
logistic and other cases, with regards to the band phases, the mean-field dynam- 
ics, and the complete desynchronization. Indeed these features are commonly 
observed in a GCM with an always expanding piece- wise linear map (i.e., the 
absolute value of the slope is always larger than 1). As an example, we plot the 
overlaid time series of the GCM (1) with 

f{x) = ax{x < 0.5); a{x — 0.5)(x > 0.5), 

in Fig. 11. Again we see single-band, 2-band, (and a higher band), and com- 
pletely desynchronized states. The dynamics belongs to the same class as the 
globally coupled tent map. 

Quite recently Pikovsky and Kurths reported on the collective dynamics of 
globally coupled tent maps with a "multiplicative coupling" |T^. They have pre- 
sented a beautiful example of clear quasiperiodic mean-field dynamics. Although 
these behaviors are somewhat similar to our globally coupled tent map, thus far 
band-type motion has not been reported in their example. 

Recall that in globally coupled logistic maps ( and in other systems with ex- 
pansion and contraction), neither clear mean- field dynamics nor complete desyn- 
chronization is found. Thus we conclude non-trivial violation of the law of large 
numbers. In the globally coupled complex Ginzburg Landau equation, collective 
behavior has some features in common with the single-band state in our globally 
coupled tent map, although it has some differences both from the logistic and tent 
maps, as for the number of clusters, and the mean-field dynamics. Classification 
of globally coupled dynamical systems waits for further studies. 

Another remaining question is about the topology of couplings. In locally 
coupled dynamical systems, collective quasiperiodic motion has been extensively 
studied p8[. The collective order of the Ising type has explicitly been constructed 
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|39| , while the antiferro type order (with a zigzag pattern) is often found in CML 
|3l| . In the former case, the collective dynamics is a fixed point, while a periodic 
or quasiperiodic collective motion is observed in the latter case. More generally, 
such quasiperiodic collective motion coexists with local chaos in the pattern se- 



lection regime of the CML |31|, as well as in the model of the Benard convection 
with a large aspect ratio [^. The collective order in the two-band state here 
can be regarded as the mean-field version of such quasiperiodic collective order 
in the zigzag regime of CML. In a CML ( with local interaction), a self-consistent 
Perron Frobenius approach has been developed |^l|. It will be important to study 
the collective dynamics in local CML by using this PF approach. Last, we note 
that the clustering in a hypercubic lattice is discussed, which lies between global 
and local interactions |^. Further studies will be necessary for the connection 



between locally and globally coupled chaotic systems. 
Fig. 1 1 

Acknowledgements 

The author would like to thank Y. Takahashi, N. Nakagawa, T. Konishi, and 
T. Ikegami for useful discussions. He would also like to thank F. Willeboordse 
for critical reading of the manuscript and illuminating comments. Thanks are 
also extended to Dr. S.Ershov for sending me his unpublished notes. This work 
is partially supported by Grant-in-Aids for Scientific Research from the Ministry 
of Education, Science, and Culture of Japan. 

References 

[1] The present article presents only a part of the talk presented at Como, 



while other parts are given in the succeeding references, as well as |22| for 



the application to biological networks. After the presentation of the talk. 



some related papers (e.g., [|T^]) have been published which are also discussed 
here. 

[2] K. Kaneko, Phys. Rev. Lett. 63 (1989) 219; Physica 41D (1990) 38 

[3] K. Kaneko, Phys. Rev. Lett. 65 (1990) 1391 and errata; Physica 55D (1992) 
368 

[4] K. Kaneko, Physica 54 D (1991) 5 

[5] for chaotic itinerancy see also: K. Ikeda, K. Matsumoto, and K. Otsuka, 
Prog. Theor. Phys. Suppl. 99 (1989) 295; I. Tsuda, "Chaotic Neural Net- 
works and thesaurus" , in Neuroconiputers and Attention, (eds. A. V. Holden 
and V. I. Kryukov, Manchester Univ. Press, 1990) 

[6] K. Kaneko, J. Phys. A 24 (1991) 2107 



11 



[7] K. Wiesenfeld and P.Hadley, Phys. Rev. Lett. 62 (1989) 1335 

[8] G. Perez et al., Phys.Rev. A 45 (1992) 5469; G. Perez, H.A. Cerdeira, Phys. 
Rev. A 46 (1992) 7492; G. Perez et al., Physica 63 D (1993) 341; 

[9] S.Sinha et al., Phys Rev. A 46 (1992) 3193; 46 (1992) 6242 

[10] D.Dominguez and H.A. Cerdeira, Phys. Rev. Lett. 71 (1993) 3359; preprint 

[11] N. Nakagawa and Y. Kuramoto, Prog. Theor. Phys. 89 (1993) 313, and 
Physica D in press; 

[12] N. Nakagawa and Y. Kuramoto, "Singularity of the Lyapunov Spectrum in 
Globally Coupled Oscillators", Physica D, submitted 

[13] V. Hakim and W.J. Rappel, Phys. Rev. A 46 (1992) 7347 and 
preprint (1994) "Noise induced periodic behavior ub the globally-coupled 
complex Ginzburg-Landau equation" 

[14] A.S.Pikovsky and J. Kurths, Phys. Rev. Lett. 72 (1994) 1644 

[15] K. Kaneko, " Information Cascade with Marginal Stability in Network of 
Chaotic Elements", Physica D (1994), in press 

[16] See also, M. Grinasty and V. Hakim, Phys. Rev. E 49 (1994) 2661; J. 
Theiler, Los Alamos Report No. LALP-92-006 (1992); W. Just, Physica D 
submitted; T. Shinbrot, preprint; T. Posche and J.A.C. Gallas, "Synchro- 
nization effects in the dynamical behavior of elevators", preprint. 

[17] K.Okuda, Physica D63 (1993) 424. 

[18] S. H. Strogatz and R.E. Mirollo, Phys. Rev. E 47 (1993) 220 S. Watanabe 
and S. Strogatz, Phys. Rev. Lett. 70 (1993) 2391. 

[19] T. Konishi and K. Kaneko, J. Phys. A, 25 (1992) 6283; K. Kaneko and T. 
Konishi, Physica D 71 (1994) 146 

[20] D.Golomb, D. Hansel ,B.Shraiman and H.Sompolinsky, Phys.Rev. A45 

(1992) 3516. 

[21] T. Ikegami and K. Kaneko, Chaos 2 (1992) 397 

[22] K. Kaneko, " Relevance of Clustering to Biological Networks" , Physica D 
(1994), in press 

[23] Nozawa, H., "A neural network model as a globally coupled map and ap- 
plications based on chaos". Chaos 2: 377 -386 

[24] A. Winfree, The Geometry of Biological Time, (Springer, 1980) 



12 



[25] Y.Kuramoto, in International Symposium on Mathematical Problem s in 
Theoretical Physics, H.Araki ed., Lecture Notes in Physics 39 (Sprin ger, 
New York 1975), p.420. 

[26] Y.Kuramoto, Physica D50 (1991) 15. 

[27] Y.Kuramoto, Chemical Oscillations, Waves, and Turbulence (Springer, 
Berlin 1984). 

[28] P.C. Matthews, R.E. Mirollo, and S. H. Strogatz, Physica D 52 (1991) 293; 

H. Daido, Phys. Rev. Lett. 61 (1988) 231 



[30] K. Wiesenfeld et al, Phys. Rev. Lett. 65 (1990) 1749; H. G. Winful and 
L. Rahman, Phys. Rev. Lett. 65 (1990) 1575; C. Bracikowski and R. Roy, 
Chaos 1 (1991) 49; K. Otsuka, Phys. Rev. Lett. 67 (1991) 1090; W.J. 
Rappel, Phys. Rev. E. 49 (1994) 2750 

[31] K. Kaneko,Prog. Theor. Phys. 72 (1984) 480; Physica 34D (1989);in For- 
mation, Dynamics, and Statistics of Patterns ed. K. Kawasaki, A. Onuki, 
and M. Suzuki, World. Sci. 1990; Theory and Applications of Coupled Map 
Lattices (K. Kaneko ed., Wiley 1993). 

[32] S. Ershov, private communication ( note unpublished, 1993, Dec.) 

[33] R. Bowen and D. Ruelle, Inventiones Math. 29 (1975) 181; D. Ruelle, Ther- 
modynamic Formalism, Addison Wesley, MA (1978) 

[34] Y. Oono and. Y. Takahashi, Prog. Theor. Phys. 63 (1980) 1804 

[35] P. Manneville, Lecture Notes in Physics 230 (1985) 319 

[36] K. Kaneko, Physica 23D (1986) 436 

[37] P. Grassberger and I. Procaccia, Physica D 9 (1983) 189 

[38] H. Chate and P. Manneville, Prog. Theor. Phys. 87 (1992) 1 

[39] H. Sakaguchi, Prog. Theor. Phys. 86 (1991) 303 

[40] T. Yanagita and K. Kaneko, preprint 

[41] K. Kaneko, Phys. Lett. 139A (1989) 47; J. M. Houlerik, I. Webman, and 
M.H. Jensen, Phys. Rev. 41A (1990) 4210 K. Kaneko, Prog. Theor. Phys. 
Suppl. 99 (1989) 263 



13 



Figure Captions 

Fig.l 

Relation between /ij„ and hout obtained from the PF equation for a single 
logistic map. a = 1.9 and e = 0.1. For reference the line /ij„ = hout is drawn. 

a) for .1 < hin < .5 with the increment 0.002. b) Blowup around the region 
with hin = hout, i-e., .34 < hin < -355 with the increment 0.00005. 

Fig. 2 Relation between hin and hout obtained from the PF equation for a single 
tent map with a = 1.9 and e = 0.1. Only the region with .585 < hin < -605 is 
plotted, while the linear dependence (with a small slope) extends to the interval 
[0,1]. For reference the line hin = hout is drawn. 

Fig. 3 Dynamics of the probability distributions Pn{x) for the PF equation with 
noise, a = 1.9, and e = 0.1. The noise is taken from a homogeneous distribu- 
tion over [— o", o"]. Successive snapshot distributions for n = 500,501,502,503 are 
overlaid, a) a = 0.000625 h) a = 0.00125. 

Fig. 4 Return Map: hn+i vs /i„ is plotted over 10000 time steps after transients. 
a = 1.5 and e = .3, obtained from a direct simulation of globally coupled tent 
map with A^ = 50000. 

Fig. 5 Lyapunov spectra for the globally coupled tent map with e = .3, and 
N = 100, obtained from the average over 5000 steps after transients, a) a = 1.5, 
a = 1.55, ■ ■ ■, a = 1.95, from bottom to top. b) The first 25 exponents for 
N = 100 and 400 (a = 1.5 ). 

Fig. 6 Return Map hn+i vs /i„ over 50000 steps (after transients), obtained 
from the direct simulation of the globally coupled tent map with A^ = 50000. 
a = 1.5 and e = .3. 

a) a = 1.5, e = .15 b) a = 1.8, e = .3 c) a = 1.8, e = .3 (by using a different 
initial condition from b)) d) a = 1.95 e = .3, A^ = 200000 , only the part with 
hn < 0.5 is plotted, while the other part for hn > 0.5 is nothing but a hne of 
hn+i = f{hn) as in c). 

Fig. 7 Return Map hn+i vs /i„ over 50000 steps (after transients), obtained 
from the direct simulation of the globally coupled tent map with N = 50000, 
a = 1.5 and e = .3. 5 examples from different initial conditions are overlaid. The 
number of elements in one group is 50000(E), 38072(D), 34563(C), 28500(B), 
25039(A), respectively. 

Fig. 8 Return Map hn+i vs hn, over 5000 steps (after transients) is plotted, 
obtained from the PF equation, for a = 1.5 and e = .3. The PF equation is 
numerically integrated with the mesh 90000. 3 examples (A,B,C) from different 
initial distributions are shown. 

Fig. 9 

Lyapunov spectra for the globally coupled tent map with e = .3, and A^ = 100, 
obtained from the average over 5000 steps after transients. 5 examples with 
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different partition to groups are overlaid. 

a) a = 1.7 h) a = 1.8: Note tlie single-band state has a slightly larger maximal 
exponent. 

Fig. 10 The (maximal) Lyapunov exponent for the PF equation, computed from 
90000 meshes, by using the average over 3000 steps. 

a)a = 1.5, while e is incremented by 0.025. Results from 3 initial distribu- 
tions. For e < .1, the distribution approaches a fixed point where the Lyapunov 
exponent is negative, while quasiperiodic behavior of the mean-field is observed 
around .1 < e < .16. 

b)e = 0.3, while a is incremented by 0.025. Results from 5 initial distributions. 
For a > 1.85, quasiperiodic behavior of the mean- field is seen. 

Fig. 11 : Overlaid time series of x„(i) of the globally coupled piecewise-linear 
map, with N = 100, a = 1.9 and e = .45, over the time steps from 10000 to 
10050. 
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